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We study dynamical heterogeneity and growing dynamical lengthscales in two kinetically con- 
strained models, namely, the one- and two- vacancy assisted triangular lattice gases. One of the 
models is a strong glassformer and the other is a fragile glassformer. Both exhibit heterogeneous 
dynamics with broadly distributed timescales as seen in the distribution of persistence times. We 
show that the Stokes-Einstein relation is violated, to a greater degree in the fragile glassformer, 
and show how this violation is related to dynamic heterogeneity. We extract dynamical lengthscales 
from structure factors of mobile particles and show, quantitatively, the growth of this lengthscale 
as density increases. We comment on how the scaling of lengths and times in these models relates 
to that in facilitated spin models of glasses. 



PACS numbers: 05.20.Jj, 05.70.Jk, 64.70.Pf 
I. INTRODUCTION 



The dramatic dynamical slowdown accompanying the 
formation of a glass is a remarkable phenomenon 0,00|. 
One explanation for the underlying microscopic cause of 
this slowdown relies on the presence of local steric con- 
straints on the movement of particles which make them- 
selves felt to an increasing degree as the temperature is 
lowered (or the concentration of particles is increased). 
Kinetically constrained lattice gas models IH 0] are 
simple caricatures of glassformers which employ local 
steric constraints as their sole means to glassiness in the 
absence of any non-trivial static correlations between par- 
ticles (for alternative thermodynamic views of the glass 
transition see e.g. 0,0)- These constrained models have 
been extensively studied (see e.g. 0,0, El El EI)- Our 
purpose here is to extend these studies to focus on the 
idea of dynamical heterogeneity 0, El El as a mani- 
festation of excitation lines in space-time and to attempt 
identification of scaling and universality classes in the 
dynamics of these models [la El El [la, El • 



The paper is organized as follows. Section ITT1 describes 
the two models we use as well as details of the computer 
simulations used to study them. Section lTTll looks at het- 
erogeneous dynamics in our models via the distribution of 
persistence times. Section llVl presents the scaling of the 
structural relaxation time and diffusion constant, the im- 
plications of which lead to a discussion of the breakdown 
of the Stokes-Einstein relation in section Section IVT1 
discusses the emergence of a dynamical lengthscale and 
gives a quantitative characterization of this length by an- 
alyzing structure factors of mobile particles. Finally, we 
end with a discussion of our results in Section IVIII 



II. MODELS AND COMPUTATIONAL DETAILS 

We present results for two kinetically constrained tri- 
angular lattice gas (TLG) models introduced by Jackie 
and Kronig |5j- These two-dimensional models are vari- 
ants of lattice models proposed by Kob and Andersen 0] . 
Each site of the triangular lattice has six nearest neigh- 
bor sites and can hold at most one particle. A particle 
at site i is allowed to move to a nearest neighbor site, i', 
if (i) i' is not occupied and (ii) the two mutual nearest 
neighbor sites of i and i' are also empty. These rules coin- 
cide with a physical interpretation of steric constraints on 
the movement of hard core particles in a dense fluid . 
We call the model with these rules the (2)-TLG because 
both mutual nearest neighbors need to be empty in order 
to facilitate movement. We also present results for the 

(1) -TLG where the constraints are more relaxed: move- 
ment is allowed as long as either of the mutual nearest 
neighbors is empty. As with other kinetically constrained 
lattice gas models, the TLG has no static interactions 
between particles other than those that prohibit multi- 
ple occupancy of a single lattice site. Therefore, initial 
configurations can be generated by random occupation of 
empty lattice sites by particles until the desired density 
is reached. 

In the computer simulations, we investigated particle 
densities, p, between 0.01 and 0.80 for the (2)-TLG and 
between 0.01 and 0.996 for the (l)-TLG. The density p = 
1 corresponds to the completely full lattice in both cases. 
For the (2)-TLG, we used a lattice with edge length L = 
128 for all densities. There exists the possibility in the 

(2) -TLG of initial configurations containing an unmove- 
able structure which percolates throughout the system 
called a backbone ]{| . Since the dynamics obey de- 
tailed balance, these backbones could never be destroyed 
in the course of the simulation. For the densities studied 
here, however, L = 128 is sufficiently large such that the 
probability of having such a configuration is vanishingly 
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FIG. 1: Distribution of local persistence times for (a) the 
(2)-TLG model (from left to right) p = 0.40, 0.50, 0.65, 0.70, 
0.77 and 0.79; and (b) the (l)-TLG model (from left to right) 
p = 0.20, 0.60, 0.70, 0.75, 0.80, 0.85, 0.92, 0.95 and 0.97. The 
abscissae are given on a scale of logarithm base 10. 



small (see @). For the (l)-TLG, we used L = 128 to 
L = 2048. A pair of vacancies in the (l)-TLG can al- 
ways diffuse freely [B| so this version of the TLG does 
not suffer from the problem of backbones in the same 
way as the (2)-TLG. At higher densities, however, one 
still needs to ensure that there are a sufficient number 
of potentially mobile particles in the system such that 
the typical dynamics of the model are observed. For the 
(l)-TLG, the number of potentially mobile particles at 
higher densities is approximately equivalent to the num- 
ber of vacancy pairs and therefore can be estimated as 
(1 — p) 2 L 2 . The system sizes at the various densities 
for the (l)-TLG simulations were chosen such that the 
number of potentially mobile particles estimated in this 
way was always approximately 100. For both models, 
periodic boundary conditions were used. 

At each density, several hundred independent trajec- 
tories of lengths 10-100 times T a , where r a is the time 
for the self-intermediate scattering function at q = (tt, 0) 
[2fjj to reach 1/e of its initial value (see below), were run. 
Trajectories were stored logarithmically for later analy- 
sis (i.e. configurations were saved after 1, 2, 4, 8, 16, 32, 
etc. sweeps). At each state point, between 128 and 256 
independent trajectories were acquired. Time was mea- 
sured in Monte Carlo sweeps. During each sweep, par- 
ticles were chosen randomly and a move was attempted. 
For the higher density runs in both models, a continu- 



ous time algorithm was used for greater efficiency |21j . 
This algorithm involved making and updating a list of 
only those particles which have the possibility of moving 
and choosing from among those exclusively during every 
move. The total time (in units of Monte Carlo sweeps) 
was then updated accordingly by adding to it the inverse 
of the number of mobile particles available during that 
continuous time step. The continuous time algorithm re- 
sulted in a speed up of our simulations by as much as 1-2 
orders of magnitude for the highest density runs in both 
the (l)-TLG and the (2)-TLG. Finally, for the distribu- 
tion of site persistence times (see below) , statistics were 
gathered over runs of very large systems (L = 1024 and 
2048). 



III. HETEROGENEOUS DYNAMICS AND THE 
DISTRIBUTION OF PERSISTENCE TIMES 

A central phenomenon behind our perspective of 
glasses is dynamical heterogeneity ^2 113 • A direct mea- 
sure of heterogeneous dynamics in glassy systems is the 
idea of persistence times |l6j,|ljj. Fig.^a) and Fig.^b) 
show the distribution of site persistence times, 7r(logr), 
in the (2)-TLG and (l)-TLG at various densities; that 
is, the distribution of times, given an initial configura- 




FIG. 2: Self-intermediate scattering function at wavevector 
q = (vr,0) for (a) the (2)-TLG (from left to right) p = 0.01, 
0.05, 0.10, 0.20, 0.30, 0.40, 0.50, 0.60, 0.65, 0.70, 0.75, 0.77, 
0.79, 0.80; and (b) the (l)-TLG (from left to right) p = 0.20, 
0.30, 0.40, 0.50, 0.60, 0.70, 0.75, 0.80, 0.85, 0.88, 0.90, 0.92, 
0.95, 0.97, 0.98, 0.99, 0.993, 0.995, 0.996. The unit of length 
is the lattice spacing. 
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FIG. 3: (a) Structural relaxation time as a function of temperature. Circles and squares correspond to the (2)-TLG and 
(l)-TLG, respectively, throughout the paper. The inset shows the same data for the (2)-TLG plotted versus T — T c and a 
power law fit (solid line) to a portion of the data, (b) Structural relaxation as a function of scaled reciprocal temperature. 
Here T g is such that r a (T g ) — 10 7 . The lines are fits to the relaxation times at low temperatures. For the (l)-TLG we use 
the Arrhenius form lnr a oc A/T, with A ~ 2.29. For the (2)-TLG we use a double exponential lnr a tx exp(a/T), with 
aw 1.76. 
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FIG. 4: The mean-squared displacement: (a) the (2)-TLG, 
p = 0.01 to 0.80; (b) the (l)-TLG, p = 0.01 to 0.996. 



tion of the lattice, before the first change at a particu- 
lar site occurs, either due to an empty site being filled 
or a filled site becoming empty. These distributions are 
multi-point functions because they depend not only on 
the state of a lattice site at the initial time and the time 
when it changes, but also on all intervening points in 
time. As has been shown in spin- facilitated models |16| . 
three distinct dynamical regimes are observed: (i) at low 
densities, there is a single peak at small relaxation times 



indicating homogeneous fast dynamics; (ii) at interme- 
diate densities, two peaks develop and co-exist, one at 
faster times and the other at slower times, indicating 
heterogeneous, fluctuation dominated, dynamics; (iii) as 
density is increased even higher, the peak at faster times 
becomes suppressed relative to the peak at slower times 
and the dynamics again become homogeneous and slow. 
In region (ii), or the crossover region, the dynamics are 
broadly distributed over several orders of magnitude in 
time. Following ^(|, we can define, qualitatively, an 
onset density, p 0l for both models where the dynamics 
begin to feel the influence of dynamical heterogeneity 
and thereby lose their mean-field character, as well as a 
crossover density p c where slow processes begin to domi- 
nate. For the (2)-TLG, p = 0.50 and (we anticipate) p c 
= 0.80 and for the (l)-TLG, p Q = 0.60 and p c = 0.85. 

In the next section, we turn to two two-point func- 
tions which have been the more conventional measures 
of glassy dynamics: the self-intermediate scattering func- 
tion and the mean-squared displacement. 



IV. DYNAMICAL SLOWDOWN 

A basic ingredient of a glassformer is a precipitous dy- 
namical slowdown over a narrow range of temperatures or 
densities. This characteristic can already be seen qualita- 
tively in the distribution of persistence times by looking 
at the movement of the mean of the distributions as p 
increases. For example, in the (2)-TLG, from p = 0.70 
to p = 0.77, the mean persistence time increase 2-3 or- 
ders of magnitude. The traditional measure of this slow- 
down is the self- intermediate scattering function, F s (q, t) 
— ( e "i'(ri(t)-ri(o))^ particularly its decay at wavevector 
q = (tt,0), Fig. |2Ia) and Fig.^b). Here, r,(t) denotes 
the position of particle i at time t. The angled brackets, 
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FIG. 5: Self-diffusion constant as a function of inverse tem- 
perature. 



(•••}, denote an average over different pairs of config- 
urations along a trajectory separated by a given time 
interval. The decay of the scattering function to 1/e at 
this wavevector is typically defined to be r a , or the struc- 
tural relaxation time, as it gives a sense of how density 
fluctuations relax at relatively short lengthscales. 

Experiments report the scaling of viscosity versus in- 
verse temperature 0- Therefore, for the particular case 
of a kinetic lattice gas, one would like to make connec- 
tions between structural relaxation time and viscosity, 
and density and inverse temperature. In both exper- 
iments and computer simulations [22l l23| . it has been 
shown that the structural relaxation time scales like the 
viscosity. If one imagines that the concentration of va- 
cancies, c, in the TLG models are like excitations or fluc- 
tuating regions of high energy, then it is reasonable to 
define a reduced temperature, T, such that c = (1 — p) 

— 1/T 

= e ' or: 



ln(l-p) = l/T. 



(1) 



The structural relaxation time plotted versus T and 1 /T 
is shown in Fig. Ola) and Fig. EJb), respectively. The 
inset of Fig. Ofa) shows a plot of r Q versus T — T c for 
the (2)-TLG where T c is taken to be the temperature 
at the crossover density, p c (as defined in section llll|l . 
and the black line is a power law fit to a portion of the 
data in the manner of mode coupling theory (MCT) [24[ • 
In MCT, T c is a critical temperature where time scales 
diverge. We see that the MCT fit is valid for almost 
five orders of magnitude in time even though there is no 
dynamical arrest at T c , as the relaxation time of the (2)- 
TLG diverges only at p = 1 Tpjj . A similar MCT power 
law fit (not shown) can be made for the (l)-TLG, valid 
for about 2 orders of magnitude. 

Fig- Elk) shows the relaxation time on a reduced tem- 
perature scale in the manner proposed by Angell 0- 
Here, T g is defined as the temperature at which r Q = 10 7 . 
We see that the relaxation time of the (l)-TLG is Arrhc- 
nius growing as r Q ~ c _A with A = 2.29, as p — > 1. The 
relaxation of the (2)-TLG in contrast is super- Arrhenius, 
and the low temperature data can be fit with a double 
exponential in 1/T 10]. In this context, the (l)-TLG is a 
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FIG. 6: Stokes-Einstein violation in the (a) (2)-TLG and (b) 
(l)-TLG. 
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FIG. 7: Fractional Stokes-Einstein exponent: scaling of the 
self-diffusion constant with relaxation time. The dashed lines 
are fits to the data at longer times. 



strong glassformer whereas the (2)-TLG is a fragile one. 

Interestingly, the time exponent A as 2.3 of the (1)- 
TLG is the same as that for the one-spin facilitated 
Fredrickson- Andersen (FA) model in dimension d = 2, 
obtained in renormalization group (RG) calculations and 
observed numcr ically [3. This suggests that the (1)- 
TLG may be in the universality class of the FA model, 
the prototypical kinetically constrained spin model for a 
strong glassformer [2^] . We will examine more of these 
dynamic scaling relations below. 

The mean-squared displacement, (|Arj(i)| 2 ) = (|rj(<) — 
rj(0)| 2 ), is shown in Fig. 0J The self-diffusion coefficient, 
D s , is defined as D s = lim t _ >00 (|Arj(t)| 2 ) / 'At. We see 
that at low densities, D s for the (l)-TLG and (2)-TLG 
coincide. At higher densities, D s for the (l)-TLG is Ar- 
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FIG. 8: (Top) Self-intermediate scattering function for the (2)-TLG at various wavevectors: (a) p = 0.20, (b) p = 0.77. For 
both graphs, from left to right: q=(w,0), (tt/2,0), (tt/4,0), (tt/8,0), (tt/16,0), (tt/32,0), (tt/64,0). (Bottom) D s (q) = l/r{q)q 2 
as a function of q at various densities: (c) (2)-TLG, (from top to bottom) p — 0.20 to 0.80; (d) (l)-TLG (from top to bottom), 
p = 0.40 to 0.996. The higher density curves for the (l)-TLG include smaller q values because of larger system sizes (see section 
II). 



rhenius, and scales as D s ~ c 2 as p — > 1, see Fig.0 This 
result is in agreement with the analytical prediction of 
[To) . Moreover, this is also the scaling of the diffusion 
constant for a probe molecule coupled to the FA model 
in any dimension [HI, further evidence that the (l)-TLG 
is in the FA model universality class. On the other hand, 
D s for the (2)-TLG is super- Arrhenius (see Fig.|SJ. The 
behavior of D a is similar, qualitatively, to that of r Q . The 
quantitative difference in their scaling with density, how- 
ever, is significant, and is an indication that relaxation 
behaviors at short and long lengthscales are not the same. 
We turn to this issue now in both the (l)-TLG and the 
(2)-TLG. 



the TLG models. We see from Fig. that the SE re- 
lation is indeed violated for both the (2)-TLG and the 
(l)-TLG, the effect being more pronounced in the fragile 
case. Moreover, we see that the densities at which the 
product D s r a begins deviating from constancy coincide 
with the onset densities, p Q 's, extracted from the distri- 
bution of persistence times. This observation reinforces 
the idea that it is the fluctuation dominated nature of 
the dynamics that leads to the SE breakdown [ri|. 

SE violation implies that the self-diffusion constant 
does not scale with the structural relaxation time as 
t^ 1 . One possibility is that it obeys a fractional SE 



V. THE BREAKDOWN OF THE 
STOKES-EINSTEIN RELATION 

An important ramification of broadly distributed het- 
erogeneous dynamics is the breakdown of mean-field 
dynamical relations such as the much studied Stokes- 
Einstein (SE) relation. This relation says that diffu- 
sion scales inversely with the relaxation time, D s r a ~ 
constant. It is a quantitative statement of the expecta- 
tion that the dynamical behavior of normal liquids should 
be similar at all but the smallest lengthscales. In super- 
cooled liquids this simple mean-field approximation fails 
[23,[13, an d' gi yen the results of previous sections, we 
would expect a similar violation of the SE relation in 




qt\p) 

FIG. 9: Data from Fig. |S| collapsed onto a master curve. 
The filled symbols correspond to the (2)-TLG and the open 
symbols correspond to the (l)-TLG. The straight line is q 2 . 



FIG. 10: Growth of mobile particle regions as a function of 
observation time At at p — 0.77 in the (2)-TLG. Black and 
grey regions indicate the location of particles and white re- 
gions indicate empty lattice sites. Particles colored in black 
have moved at least one lattice spacing in a time At whereas 
particles colored in gray have not. (Top, from left to right: 
At = 10 3 and 10 4 ; bottom, from left to right: At = 10 5 and 
10 6 ; r a ~ 10 5 at this density). 



law, D s ~ where £ < 1. This is observed in ex- 
periments [23 , and is obtained theoretically for probe 
diffusion in the FA and East models [13 (see also (29|). 
Figure shows that the diffusion constant also obeys a 
fractional SE law in the TLG models. The SE expo- 
nent is £ w 0.88 for the (l)-TLG, which is the value ex- 
pected for the FA model in d = 2, £ « 2/2.3 [l7|. In the 
case of the (2)-TLG, despite the fact that D s and r Q are 
both super- Arrhenius, we find that the scaling exponent 
is temperature independent at large densities, £ w 0.58. 
The deviation of this exponent from 1 is larger than that 
for both the FA and East models in two dimensions ^3] • 
It indicates a larger violation of the SE law, consistent 
with the fact that the (2)-TLG is more fragile than either 
of those models [Toj . 

VI. DYNAMICAL LENGTHSCALES 

A. Indication of a dynamical lengthscale from a 
two-point function 

Since the growth in timescales and the violation of the 
Stokes-Einstein relation in the TLG models are clearly 
not tied to a growth in static lengthscales, we turn now to 
the discussion of dynamical lengthscales. Such a length- 
scale can be inferred from examining the relaxation be- 
havior of the self-intermediate scattering function over 
more than one wavevector at different densities. One can 
appreciate this fact qualitatively by looking at the decay 
of F s (q, t) for the (2)-TLG over several values of q at low 
and high density, Fig. |§J (top), [l^. At low density, the 



FIG. 11: Same as Fig. EH for the (l)-TLG. Here, p = 0.95. 
(Top, from left to right: At — 10 and 10 2 ; bottom, from left 
to right: At = 10 3 and 10 4 ; r a ~ 10 3 at this density). 



decay of the various curves looks similar at all wavevec- 
tors (except for the largest wavevector) whereas at high 
density, even the curves at intermediate wavevector differ 
greatly from the simple exponential form seen at smaller 
wavevectors. The high density curves bunch up at inter- 
mediate to large q indicating that the relaxation behav- 
ior at these lengthscales is different (i.e. slower) than 
one would expect from the behavior at larger length- 
scales. Similar behavior has also been observed in the 
Kob-Andersen kinetic latti ce g as model and kinetically 
constrained spin models 0, Il9| . 

To quantify the above behavior, we proceed as 
in |3(i l3l| . In the hydrodynamic regime, we have 
lim 9 ^o Fs (<3S t) ~ exp(— D s q 2 t), and one expects the 
product D s r(q)q 2 to be independent of q, where r(q) 
is the time when the intermediate scattering function 
at wavevector q decays to 1/e. In Fig. |S] (bottom) , we 
plot the quantity D s (q) = l/r(g)g 2 as a function of q at 
various densities. A flat line independent of q indicates 
normal diffusive behavior whereas a downward bend sig- 
nifies a change to sub-diffusive behavior. As density in- 
creases, the curves begin to bend at a smaller and smaller 
wavevector, q* . This behavior is indicative of a gro wing 
dynamical lengthscale as density is increased [l9l l30l l3ll| . 

Following the prescription of , we extract a length- 
scale, I*, from Fig. El as t ~ y/D s (q -> 0)r Q . This 
lengthscale determines the onset of Fickian diffusion. Us- 
ing i* to rescale space, and using r{q)q 2 to rescale time, 
the data from Fig.[S]can be collapsed onto a master curve, 
Fig. El 



7 






FIG. 12: (Top) Growth of dynamical heterogeneity length as a function of observation time, At, for (a) the (2)-TLG, p = 0.70 
to 0.80; and (b) the (l)-TLG, p = 0.95 to 0.996. (Bottom) (c) Scaling of the maximum value of the dynamical heterogeneity 
length, Imax, from (a) and (b) with relaxation time. The inset shows l m ax versus c for the (l)-TLG, on a log-log scale to 
highlight power law scaling. 



B. Direct observation and quantification of a 
dynamical heterogeneity lengthscale 

We can study dynamical lengthscales in the TLG mod- 
els directly by observing a trajectory over a time At and 
coloring particles which have moved at least one lattice 
spacing. Snapshots of applying this procedure to trajec- 
tories of the (2)-TLG and the (l)-TLG at high particle 
densities over progressively longer At's are shown in Fig. 
IHfland Fig.m Mobility is indeed correlated: mobile par- 
ticles are clustered and the clusters of mobility at earlier 
At act as seed particles from which subsequent mobility 
grows. Moreover, there is a qualitative difference in the 
shape of the clusters in the (2)-TLG and the (l)-TLG. 
In the fragile model, the clusters are smooth and more 
anisotropic, indicating a directed growth of mobile re- 
gions. In the strong model, the clusters are rugged and 
isotropic. These observations can be understood as aris- 
ing from the difference in the local constraints of both 
models. The strict two site facilitation rule of the (2)- 
TLG requires cooperative, hierarchical rearrangement of 
particles for movement whereas the one site facilitation 
rule of the (1Y-TLG allows for the random diffusion of va- 
cancy pairs |fj . These same correlations between fragility 
and the smoothness of interfaces and between slow and 
fast dynamically heterogeneous regions are present in 
other facilitated models |l l 111 Ip ] . 

To quantify the above ideas, we extract a dynamical 



lengthscale, l(At), from structure factors of the mobile 
particles 0, 0, I^O, Hl|. Motivated by the mobility 
criterion described in the previous paragraph and de- 
picted in Fig. ^]| and Fig. ^] we consider the binary 
field n r (t; At) = p r {t)[l — p T (t + At)}, which gives a sig- 
nal only when there is particle motion at r over the range 
of time between t and t + At. Its structure factor is a 
four point function: it measures a correlation function 
which depends on two points in time, t and t + At, and 
two points in space, r and r'. 

We define the structure factor for the mobility field as 
the following normalized correlation function: 

1 (^A^^ 



L 2 



(<5n r (t;At) 2 } 



where Sn r (t; At) = n r (t; At) — (n r (t; At)) is the deviation 
of n r from its average value and <5n q (t; At) is the Fourier 
transform of 5n T (t; At): 

<5n q (t; At) = ex P ( T? r ' q ) <M*! At). (3) 



The angled brackets, (•••), denote an average over differ- 
ent pairs of configurations along a trajectory separated 
by a given time interval At. We then define the length- 
scale, (-{At), to be proportional to the inverse of the first 
moment, qAt, of the circularly averaged structure factor, 
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S(q n ;At) 37]. That is, i(At) ~ l/q At where: 

qAt = J2^S(q n ;At)/J2S(q n ;At), (4) 

n n 

Here, g„ = 2im/L and n = 0, 1, 2, i/2. The 
lengthscale was normalized such that t extracted from 
the structure factor of a random configuration of parti- 
cles on the lattice (i.e. an ideal gas) was unity. 

Fig. H2f a) and Fig. I12f b) show £(At) at various densi- 
ties for the (2)-TLG and the (l)-TLG. The basic shape 
of these curves is as expected: as At — > 0, mobility is 
sparse and uncorrelated so £ approaches unity and as 
At — » oo, everything becomes mobile and I once again 
tends towards unity. In between, as the pictures in Fig. 
ITU1 and Fig. ^] suggest, mobility clusters together and 
grows. Looking at the maximum of these different curves, 
£ max (At max ) , a growing lengthscale is clearly evident as 
p increases. It is important to note that, in general, 
£max 7^ £* 0- Fig. 0c) shows the value of £ max 
plotted versus the structural relaxation time for both the 
(2)-TLG and the (l)-TLG. At short relaxation times, the 
curves merge and approach the ideal gas value of one. As 
relaxation times increases, the dynamical heterogeneity 
lengthscale for the strong version of the model is always 
lar ger than that of the fragile version at a fixed value of 
r a [l5|. We also find that the observation time, At max , 
at which the maximum lengthscale, £ ma x, occurs, scales 
with the structural relaxation time, r Q , for both models 
(not shown). 

If the (l)-TLG is in the universality class of the FA 
model, then we would expect £ max to scale as a power 
of both the excitation concentration, £ rnax ~ c"i an d 

1/z 

of the relaxation time, £ max ~ r Q . This appears to be 
the case, as shown in Fig. I12f c). For the correlation and 
dynamic exponents we find v « 0.89 and 1/z rj 0.36, in 
reasonable agreement with Ref. [l8|, i/«0.7 and 1/z = 
v/ A w 0.3 [34(. The z exponents shown in Fig. ll2l' c') are 
what we would expect from Fig. I12f b) where a range of 
£(At) curves at different densities merge at early times 
and display power law scaling with similar exponents. 

As alluded to earlier, the dynamics of the (l)-TLG 
at high densities is controlled by the motion of vacancy 
pairs. The physics of these vacancy pairs is similar to 
excitations in the FA model. Vacancy pairs have the 
ability to interact with other lattice vacancies in order 
to branch and coalesce. It is important to note that 
evidence of these interactions can only be seen in sim- 
ulations of large enough system sizes where the number 
of vacancy pairs is approximately 50-100. As mentioned 
in section ITT1 this requirement leads to system sizes, for 
example, of L = 2048 for p = 0.995. 

Fig.ll3lshows the structure factors of the mobility field, 
equation J2|, measured at the structural relaxation time, 
At = r a for the (2)-TLG and the (l)-TLG. The curves 
are scaled in a manner suggestive of a coarsen ing process. 
The collapse of the various structure factors |3S| implies 
that, in the glassy regime, increasing density corresponds 




0.01 0.1 1 10 100 



FIG. 13: Structure factors of the mobility field, equation 
measured at At — r a , for various densities in the (a) (2)-TLG 
and the (b) (l)-TLG. The axes are scaled as indicated. 



to a coarsening of dynamic heterogeneity fields in space- 
time. An appreciation for the self-similarity of dynamic 
heterogeneity fields at different temperatures has already 
been noted in spin- facilitated models (lH fl8| . 

Both sets of scaled structure factors in Fig. ^5] have 
an intermediate power law regime going as fc~ 2 ' 3 for the 
(2)-TLG and £r 2 for the (l)-TLG. One explanation for 
the difference in exponents could be the following. Porod 
law scaling, k~( d+1 ^ , arises from a system which is exten- 
sive in interfaces. In two dimensions, this implies that a 
system with a scaling exponent closer to 3 would have 
smoother interfaces. This interpretation is consistent 
with the snapshots of the dynamic heterogeneity fields 
in Fig. [inland Fig. [HI 

The structure factors for the (l)-TLG imply a value 
of the dynamical exponent 77 very close to zero, S(k) ~ 
l/k 2 ^ 11 as k — * 00. This gives a prediction for the expo- 
nent 7 via the scaling relation 7 = (2 — r))v± of 7 w 1.8. 
The 7 exponent controls the scaling of the dynamic sus- 
ceptibility, x = S(q — > 0; p, r a ) with the concentration 
of excitations. The inset to Fig. ll3f b}shows that this 
expectation is approximately satisfied |l8j . 



VII. DISCUSSION 

Despite their simplicity, the constrained lattice gas 
models we have studied show the essential features of 
glass forming liquids, such as a precipitous dynamical 
slowdown and dynamical heterogeneity. They are inter- 
mediate between the fully coarse grained kinetically con- 
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strained spin models such as the FA and East models, and 
atomistic models such as the binary Lennard- Jones mix- 
ture. We find a broad distribution of persistence times, 
especially in the (2)-TLG (Fig. QJ. From the scaling of 
the structural relaxation time it follows that the (2)-TLG 
is a fragile model and the (l)-TLG is a strong one, consis- 
tent with the predicted scaling of the self-diffusion con- 
stant in [lOj . Fragile behavior versus non-fragile behavior 
coincides with hierarchical versus diffusive propagation of 
excitations IM|35l , and the former follows from directional 
persistence 15] as evident from the patterns of dynamic 
heterogeneity seen in Fig. 1101 and Fig.^] Dynamic het- 
erogeneity produces length-time scaling and decoupling 
phenomena. Dynamic heterogeneity is present in both 
strong and fragile materials, not only in the latter. This 
is consistent with recent molecular dynamics simulation 
on silica [3(j and earlier theoretical predictions [TR HtH . 
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